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Abstract 



The lepton and dilepton charge asymmetries from Bd and B s are predicted to be small 
in the standard model, whereas new physics could increase their values significantly. In 
this paper, we explore the use of the lepton asymmetries as a probe of the flavor structure 
of supersymmetric theories. In particular, we determine the sensitivity to parameters of 
various models. We find that in many interesting models which attempt to address the 
supersymmetric flavor problem, the mixing structure is such that it could be possible to 
detect new physics. The predictions are model dependent; with a measurement in both 
the B s and Bd systems one can hope to constrain the flavor physics model, especially once 
squarks are detected and their masses measured. Thus, lepton charge asymmetries can be 
used as an alternative means of searching for new physics and distinguishing among potential 
solutions to the flavor problem. They are interesting precisely because they are small in the 
standard model and are therefore necessarily evidence of new physics. 



1 Introduction 



The next few years will be an exciting era for B physics, with the detailed investigation of B 
hadrons at 6-factories. Particularly exciting is the potential for studying CP-violation in the 
B system, both within and beyond the standard model (SM). This affords the opportunity 
to look for new physics, and in fact might yield the first evidence for physics beyond the 
standard model. Hopefully this new physics will be further studied directly so that we will 
establish its origin. Whatever this new physics proves to be, a detection at fe-factories should 
give new information which will not be accessible to high energy colliders. For example, 
should this new physics prove to be supersymmetry, detailed studies in the B system give a 
unique opportunity to probe the flavor structure of extensions of the SM. 

Non-standard model CP-violating effects could be revealed by testing whether measure- 
ments agree with the SM allowed range. Processes for which the SM contribution vanishes 
(or is negligibly small) offer an important complement to these studies. In this case, any 
observation or non-observation of CP-violation can be interpreted directly as a constraint 
on physics beyond the SM. From this point of view, a measurement of the dilepton or of the 
lepton charge asymmetry is of great interest. The dilepton charge asymmetry is defined asQ 
1, J, |, |, § " 

where l ++ [/ ] denotes the numbers of / + / + dilepton pairs observed. In the 6-factories, 

they come from the mixing and decay of the BB pairs, while in the hadron colliders, the final 
hadronization states can be any combinations of B + , B^, B s , Aj, and their conjugates. In the 
absence of CP- violation, this quantity clearly vanishes. In early studies of the dilepton charge 
asymmetry J7], the KM angles and top quark mass were not sufficiently well determined 
to be certain that a measurement in excess of 10~ 3 would signal new physics. As we will 
see, this quantity is now determined to be small in the SM for both Bd and B s , but can be 
significantly larger in non-standard flavor models. 

Another useful quantity to look at 0, ||, ||, || is the total lepton charge asymmetry l^, 
which is defined by 

l± = FTF- ™ 

Here = N(BB — > Z +( ~)) is the total number of positively (negatively) charged primary 
leptons coming from the decay of bb pairs. This quantity is smaller than An, but should be 
measured with better statistics. 

These lepton asymmetries are sensitive to the phase difference between T 12 and M 12 . In 
the SM, the dominant contribution to these quantities has the same phase, and is therefore 
suppressed. The dominant source of enhancement in non-standard physics is a new contribu- 
tion to mass mixing, which would generally carry a different phase from the standard model 
contribution. In models for which the new source of mixing is comparable to that of the SM, 
one can expect a substantial change in the prediction of the lepton charge asymmetry. We 



1 Here we use the convention that B contains a b quark, thus decays into a l + if there is no mixing. This 
is opposite to the convention used in Q. 
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will show that there are many models for which one would expect up to an order of magni- 
tude enhancement over the SM prediction for Bd and two orders of magnitude enhancement 
for B s . In fact, because the predictions for asymmetries due to CP violation in the B s and 
Bd systems is different, one can hope to use this measurement to help distinguish among 
potential solutions to the supersymmetric flavor problem. For this reason, we take an explic- 
itly model-dependent approach to our results, and explore the predictions in various models 
already existing in the literature. They do not necessarily include the real world solutions 
to the flavor problems, but are nonetheless sufficiently general to illustrate the usefulness of 
the lepton asymmetry methods. 

There are several ways to measure the lepton asymmetries. One can measure both the 
single lepton and dilepton asymmetries at the dedicated B factories. These will of course only 
be sensitive to new physics in the B& sector. It would be extremely interesting to complement 
this measurement with the measurement of lepton asymmetries at a hadron collider, which 
will be sensitive to the asymmetries in both B s and Bd- With all such measurements (or 
even some fraction of them) one should be able to distinguish new mixing contributions to 
either the B s or Bd systems. We will see in particular that many supersymmetric flavor 
models yield a large deviation for at least one of the above systems. 

Because the SM prediction for Bd is small, and even smaller for B s (we discuss how 
small later), any lepton asymmetry measurement in excess of this value is a clear signal 
of new physics. Because of the suppression from AT /AM, a sensitivity of at least 10~ 2 is 
essential. So any asymmetry within this range will be an important signal of new physics. 
We emphasize that even without flavor tagging, a measurement of a CP-violating asymmetry 
in excess of the SM prediction will be an exciting signal of new physics, which ultimately 
complementary measurements should disentangle. 

The importance of studying lepton asymmetries has been considered previously in |1 



specifically in the context of B s , and in |K| for Bd- Re(e#) (which can be related to An by 



An ~ 4Re(e#) |TT|) in left-right symmetric models has also been discussed in ||12| . However 



no one has as yet done a detailed study of the potential significance of this measurement 
for distinguishing models for the supersymmetric soft masses or combined the information 
from both Bd and B s . Of course, any non-standard model can be studied in the light of 
the measurement of the lepton charge asymmetry and thereby constrained. For specificity, 
and because of its likelihood as the source of a non-standard model CP- violating effect, we 
chose to study the specific case of contributions from soft scalar masses in some non-standard 
models of flavor physics that exist in the literature. We find that many squark mass models 
designed to address the flavor problem in supersymmetry will give rise to a significantly 
larger lepton asymmetry in either the Bd or B s system for reasonable parameters. A simple 
order of magnitude estimate shows that the box diagrams would be comparable (assuming 
mixing similar to CKM mixing) for superpartner masses of order a few hundred GeV. 

We begin in Sec. ^| by presenting the basic formulas relevant to BB mixing. We ex- 
press the dilepton asymmetry in terms of (Ar/AM)sM> and the relative amplitude and 
phase difference between the supersymmetric and SM contributions. We then briefly review 
why supersymmetry gives new contributions to neutral 5-meson mixing through non-trivial 
squark mass matrices. In Sec. [| we give a brief review of three soft supersymmetry-breaking 
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scenarios which have been devised to address the flavor problem in supersymmetry mod- 
els: alignment, non-abelian flavor symmetry, and heavy squark models. These suggestions 
solve or relax the Flavor- Changing Neutral- Current (FCNC) and CP-violation problems. In 
Sec. £|, we first present the model-independent lower limits on the phase difference between 
Mf 2 USY and Mp 2 M and the rriq range (linear with respect to the mixing angle), assuming the 
experimental sensitivity of measuring either dilepton asymmetry is 2 x 1CT 3 . For any par- 
ticular ansatz for the squark mass matrices, this result can be interpreted as a sensitivity to 
rriq and mixing angles. We show the parameter ranges of the models of Sec. ^|, which can be 
probed with the measurement of the dilepton asymmetry. In Sec. |5|, we conclude. 



2 Dilepton Charge Asymmetry 



For the B B basis, one has the Hamiltonian 



The mass eigenstates are 

with eigenvalues 
where 



/ m — ir/2 m 12 - ir 12 /2 \ 

I M* 2 -ir* 2 /2 M-ir/2. J 1 • 



Bi,2= i = (\B)±V\B)), (2.2 
M 1>2 - ir 1>2 /2 — M — ir/2 ± A, (2.3 



V 



M{ 2 - irt 2 /2 
^M 12 -ir 12 /2' 



(2.4 



A = y(Mi2 - ir 12 /2)(M 1 * 2 - ir* 12 /2). (2.5 
The quantities r,f are defined as p|, |13| | 

_ Pb^b i + y 2 (9(] 

r ~ Pb^b H 2 2 + x 2 -y 2 ' 1 

p b ^b = I i2 x2 + y 2 (27 

Pb^b m 2 + x 2 -y 2 > 

where x = AM/T, y = AT/2T, AM = M x - M 2 = 2 Re A, r = {V x + T 2 )/2 and Ar = 
r x - T 2 = -4ImA. 

When a BB pair is produced, it can mix and later decay into Z + Z + or . Thus, we can 
replace the Z ++ and I in Eq. ( |1.1| ) by N(BB) and N(BB), which is the number of BB 

2 The r, f defined here is the interchange of r, f defined in ^ [l3| because of the opposite convention used 
in defining B. 
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(BB) final states observed in a sample of events from a process where a BB pair is initially 
produced f2j. The dilepton asymmetry can then be written as [§. [| 

N(BB)-N(BB) _r-f _ \y\ 4 - 1 _ Im(r 12 /M 12 ) ^ / T 12 \ 
w N(BB)+N(BB) r + f \r)\ A + 1 1 + l/4|r 12 /M 12 | 2 ~ VM 12 A l ' J 

The last approximation holds if |ri 2 /Mi 2 | <C 1, which is the case for the BB systems even 
in the presence of new physics [ 14] . 

This formula is true whether or not BB is produced coherently. However at a hadron 
collider when a B d is not necessarily produced in conjunction with a B d , one needs to account 
for all possible fragmentations. In this case we derive 

A _ l ++ ~ I _ feW, + ^Xsfs) A++ - (^-X d f d + ^XJ S ) A- 



where 



and 



, \ ,. \ Xdfd , X s f s 



A = X f + X f H — _i£z — | — _^>j_^ / 21Q) 

1 + r d 1 + r s 



A = X f + X f H — — — ^— — I — _±±£_. (2.H) 

1 + r d 1 + r s 

Here / l is the probability to hadronize as a state % and A* is the leptonic branching fraction. 
It is readily seen that when working to leading order in CP asymmetries, the formula (|2.9| ) 
reduces to 

(AOtotai = (A u ) d A d2 + (A u ) s A s2 , with (A u ) d , s = rd ' s ~ r ^\ ( 2 .12) 

rd,s + r± s 



where 



\ d f d r d (A + / + + A / + X d f d + X s f s ) 

d2 (1 + ^)MA + / + + AA/a + a^ + a^ )( ^ + M^)' ^ > 

A XsfsTs (A+/+ + X A f A + X d f d + X s f s ) 

32 (l+r s )2(A+/+ + A A / A + ^ + ^)(^-^^ " ' 



l+r s 



and we define A 2 to be the ratio between A s2 and A d2 , which gives 



A s2 X s f s r s (l + r d ) 
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A-2 = —r- = : bJ : Jr-^r (2-15) 
A d2 X d f d r d (l + r s ) 2 

The most recent results for f and A* are /+ = 0.39 ± 0.04 ± 0.04, f A = 0.096 ± 0.017, 
f d = 0.38 ± 0.04 ± 0.04, f s = 0.13 ± 0.03 ± 0.01 |T§, A+ = (10.3 ± 0.9)%, A A = (9.0±H)%, 
X d = (10.5 ± 0.8)% and X s = (8.1 ± 2.5)% 0]. For the BB system, \rj\ in Eqs. and 
is very close to 1, y/x = AT/2AM is very close to 0, thus 



x 2 



■-■-2 + ^ (2 - 16) 

For B s , x is large because of large mixing; thus B s is almost 100% mixed and r — > 1. For 
-B,i, x = 0.734 ± 0.035[|16], which gives r = 0.21 ± 0.02. Putting together these numbers, 
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we find that Ad2 = 0.53, A s2 = 0.25, A 2 = 0.46; this means that, although B s is 100% 
mixed, its contribution to the dilepton rate is less than B^s because the number of B^s 
produced is about three times larger than B s . In our analysis below, we will consider 
super symmetric contributions to both B^ and B s mixing. It should be borne in mind that 
a better measurement of An is required to achieve the same sensitivity to the parameters 
relevant to B s . 

The total lepton charge asymmetry has a different form when expressed in terms of r 
and f in the case of a coherently or incoherently produced BB pair. When it is produced 
coherently, for example, in T(4S f ) — > BdBd, the total lepton asymmetry is given by |@, 0, H H 

= l ^ zJ r = ( 2 - 17 ) 

/+ + /- 2 + r + f y J 

One can simplify the total lepton charge asymmetry by observing 

l± = - ^,,2,''" , :„ , : « 0.17 (A u ) d for B d , (2.18) 
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here again we take \rj\ — > 1 in the denominator. In the 6-factories, only BdBd is coherently 
produced. 

When the BB is produced incoherently, the total lepton asymmetry becomes || 

I* = - — , (2.19) 

1 + r + r + rr 

which can also be simplified as 

' ± = 'I"" « { °„ 2 ' ft* !" (2.20) 
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0.5 (A u )s for B a . 

The single lepton asymmetry measures the same quantity, Im(ri 2 /Mi 2 ), as the dilepton 
asymmetry. The prediction for the dilepton asymmetry is bigger; however because both 
leptons must be tagged, the statistics are smaller. Yamamoto || has argued that the single 
lepton asymmetry will give a better measurement at a dedicated 6-factory. At a hadron 
collider, the contribution from both B d and B s as well as other sources of single leptons 
must be accounted for: 

(/ ± )total = l^Adl + ifAsi, (2.21) 

where 

A _ ^ d f d A — ^ s f s O 99^1 

dl A+/+ + A A / A + \ d fd + Xsfs ' sl + A A / A + Xdfd + Kfs ' 1 ] 

and l^g takes the form of Eq. Q2.19| ). If the leptonic branching ratios were the same for all 
6-hadrons, Adi, s i would reduce to fd, s - We define A\ to the the ratio between A s \ and Adi, 
which is 

A x = — = ^ s ^ s (2.23) 
A d i Xdfd 
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For the values of fd, s and Xd, s given before, we get A& = 0.40, A s \ = 0.11, A% = 0.26. 
However, to measure Imr 12 /M 12 with the same sensitivity for Bd and B s requires the same 
factor as with the dilepton asymmetry, in light of Eq. ( |2.20| ) . 

In the SM, the phases in Ti 2 and Mi 2 are approximately equal. Thus, 

AM SM ~2|M 1 S 2 M |, Ar SM ~2|r^ 2 M |. (2.24) 

The SM contribution to the dilepton charge asymmetry is generally small |4], ^, [H], [17]] 

f 10~ 3 for B d , 
\M ~ | 10 _ 4 for Bg (2-25) 

The current preferred solutions of the unitary fits to the CKM matrix yield: p ~ 0.12 and 
T] ~ 0.34 JI8|] , which gives \Au\ ~ 4.2 x 10~ 4 for Bd in the SM. When the mixing angles vary 
over their currently allowed range |n|, the dilepton asymmetry can be as large as 1.9 x 10~ 3 
for Bd- Until the angles are better known, a lepton asymmetry in excess of this number is 
required in order to test new physics, so in this paper this will be our benchmark. That 
is, we assume the experimental sensitivity is good enough to measure down to the largest 
possible standard model value. However, once these angles are better determined, even 
smaller values could indicate non-standard contributions. We will show that there are many 
interesting models that predict a contribution at this level. The current experimental bound 



on the dilepton asymmetry is \Au\ < 0.18 |20, |T] at the 90% confidence level, far below the 



level of interest. We stress the importance of better measurements at Run II of the Tevatron 
and at 6-factories, with the ultimate goal of at least this sensitivity. 

One should bear in mind the reduced contribution of B s relative to B d (assuming equal 
Im(ri 2 /Mi 2 )) which means better experimental sensitivity is required to study B s at the 
level assumed. This would require decoupling any possible standard model contribution 
which could be present. 

Another caveat is that the small standard model rate is based on a quark calculation, 
and relies on a sensitive cancellation between intermediate states with light up-type quarks. 



Wolfenstein |22[ has argued that the quark model calculations might not be reliable and 
predicts a much larger rate based on 100 % violation of duality. Even with a smaller violation 
of duality, of order 20 % for the c c intermediate state, we find the standard model prediction 
could be increased by a factor of 3 if the quark model rate is an overestimate. However, 
there is no evidence as yet for this violation. It will be interesting to better test the 



assumption in the future by a better measurement of n c , the average number of charm (or 
anticharm) quarks in the hadronic final state of a B decay. 

As is well known, there can be large FCNCs in supersymmetric models because of the 
many new potentially flavor-violating parameters. In particular, the squark mass matri- 
ces introduce the possibility of new flavor- violating effects. These effects can be described 
through the mass matrices M"[ L , M^j, and M\ R . Because of potential new contributions to 
BB mixing, the phase of T 12 and M 12 should be different. Assuming that supersymmetry 
does not substantially change T 12 (since it only contributes at higher order) and defining 
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Figure 1: Plot of dependence of A u /(AT/AM) S m on h, with 9 = vr/4, vr/2, 3tt/4 



AfSUSY^j^-SM _ ^ e w ^ dilepton asymmetry A// is 

/ r?^ \ _ / Ar \ hsme 

11 ~ \mWTmW^) ~ VAmJsm l + 2/icos0 + fc 2 ' ^ 2 - 26 ) 

which is similar to a formula that was also presented in Ref. ||. In the SM, (Ar/AM)sM is 
small ||§§: 

/ AT_\ f (1.3 ± 0.2) x 1(T 2 for B d , , , 

\AJW7 SM 1 (5.6 + 2.6) x 10~ 3 for B a . { } 

The errors break down as follows: in Bd, ±0.1 coming from = 4.8 ± 0.2 GeV; ±0.1 
coming from m t = 165 ± 6 GeV; ±0.07 coming from the CKM matrix elements = 
0.003 ± 0.0008, \V cd V; b \ = 0.0086 ± 0.0007, \V td V* b \ = 0.0084 ± 0.0018 JTB1; ±0.02 coming 



from ?7qcd = 0.55 ± 0.01 1 26| . In B s , the dominant error comes from B$/B. Here Bs and 
B are the "bag" parameters used to estimate the matrix element Qs = (biSi)s-p(bjSj)s-p 
and Q = {biSi)v-A(pjSj)v-A respectively (see Ref. @). There are ±2.3 from Bs/ B varying 
between 0.7 and 1.3; ij; 7 from varying /i between m&/2 and 2m^; ±0.4 from mj and ±0.4 
from m t |] The dilepton charge asymmetry An can be enhanced over the SM value by the 



second factor in Eq. ( 2.26 ). Notice that this factor reaches its maximum when h = 1, which 
gives 

(Ar \ i 9 

This dilepton asymmetry is larger when 9 is larger, especially when 9 is close to n. We plot 
the dependence of An/ (AT / AM)sm on h by assuming the CP- violating phase difference 9 
to be 7r/4, vt/2, 3vt/4 in Fig. [I]. Notice that An is heavily suppressed if h is either very large 
or very small. This is because AM is large when h is large, while there is no significant 
new contribution to AM when h is small. For a given experimental sensitivity (A^) min , 



The mt taken in Ref. 24 is m t — 176 ± 9 GeV. It is slightly different from the mt we used in our 



calculation, which is the more recent experimental result |16 
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there are corresponding h m \ n and h max (or none if 9 is too small) to which the measurement 
is sensitive. This in turn gives (mj) min and (m 5 -) miB by the formulas given in Sec. [| for 
particular mixing angles. A measured dilepton asymmetry would constrain nig to be in a 
range between (m g -) min and (ro ? -) max . The precise numerical results will be presented in Sec. [|. 

In the literature, there are two different parametrizations of the effects of non-trivial 
squark matrices. One can diagonalize the gqq coupling and quark mass matrices while 
keeping all the mixing effects in the q propagators. This is called the "mass insertion" 



method p7[ . One can also work in the mass eigenstates of quarks and squarks with off- 
diagonal gluino couplings (we call this the "vertex mixing" method), and consider only the 
contribution from the lightest generation. The "mass insertion" works well when the squarks 
are near degenerate, that is, for m 2 = m 2 (l + Xj), i=l-6 for q^ and q^, Xi <C 1. Here fh is 
the average squark mass P8| . When the squark masses are not so degenerate but have the 
same order of magnitude, the mass insertion method can still be a good approximation if 
the average mass m is chosen appropriately. The "vertex mixing" is a better approximation 
when one generation is much lighter than the other two since the contributions from the 
heavy generations are suppressed by their masses. Notice that the mass insertion method 
assumes a GIM-like cancellation of the leading term, which is why the results using vertex 
mixing and mass insertion can be different. We will use both methods in our numerical 
calculation below, according to which is more appropriate. 

In the flavor eigenstate basis of both quarks and squarks, the mass matrices for the up 
and down sector quark and squark are M u , M d , M u2 and M d2 , where the 6x6 squark mass 
matrix M d2 can be written in terms of 3 x 3 matrices Mmn {M, N = L, R): 

""(IS)' (2 ' 29) 

The off diagonal M RL and Mf 2 R are usually very small due to the suppression by mz/fh 
and quark masses (in particular for our purposes there is a A;, suppression). In addition, 
the decay rate of b — > constrains (Slr)23 (will be defined below) to be smaller than 



1.6 x 10 (mg/500GeV) 2 [28[|. We neglect Mf^ and diagonalize the squark mass matrices 



Mf L and M RR in the 

mass basis of quarks and squarks, which defines the mixing angle V 

by 

V d M d V d+ = dmg(m d ,m s ,m b ), (2.30) 
V£M»V%* = diag(m u , m c ,m t ), (2.31) 

ViMf L Vt = diag(m 2 L ,^L^L), (2-32) 
V^M d2 R V d+ = dmg(m 2 dR ,m 2 sR ,m 2 bR ). (2.33) 

The gqq vertices are in general not diagonal: the coupling mixing matrices (analogous to the 
standard CKM matrix) are 

K d = V d V d+ , K d R = V d V d+ , (2.34) 

and similarly for the up system. This method of calculating flavor- changing effects is par- 
ticularly useful when the mass eigenstates are very non-degenerate. 
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We can also work in the basis where the gqq couplings and quark mass matrices are 
diagonal. All the mixing is now in the squark propagators, which can be expressed in terms 
of the dimensionless parameters (S^mn 

c _(Sll Slr\ J_(V£MtlV£ + VgM&V* 
MN \6rl S rr ) m*\vZM&Vt V«Mi R V% 

The mass insertion method is valid if (S^mn (i ^ j) is small. 



(2.35) 



3 Models of Soft Super symmetry Breaking 

As is well known, the many additional parameters in supersymmetric models can introduce 
large dangerous FCNC and CP-violating effects. The parameters must be such that the 
experimental bounds on e^, K, B, D mixing, the electric dipole moment (EDM) of electron 
(d e ) and the neutron (d n ) and branching ratio of b — > 57 |28], |3(| are preserved. Different 
scenarios have been proposed to solve, or at least relax the FCNC and CP-violation problems 
3jJ . It is possible that none of these are the true solution, but they serve as useful straw 
men. In this paper, we will discuss three of them: alignment |52|, |33|, |34], non-abelian 
0, ||, ||, |g, H, H, H and heavy squark models @, ||, ||. This list of 



models ||, g§, g7|, 

references is incomplete but incorporates the models we study. Any model can be interpreted 
as we do with these. Models with nearly exact universality, such as gauge-mediated models, 
are of course an intriguing possibility for solving FCNC problems; however, FCNC effects are 
generally suppressed, and we therefore do not discuss these here. In particular, the phase of 
T 12 and M12 would still be correlated if universality were assumed as a boundary condition. 
Should the effects we describe be observed, gauge-mediated models would be excluded. 



The idea of alignment |32|, |33|, |34| is that the squark mass matrices are aligned with the 
quark ones so that the Kl jR in Eq. ( p.34|) are close to the identity; that is, the off-diagonal 
terms are small. Therefore, there are no large contributions to FCNC. Such models can be 
constructed with an abelian horizontal symmetry H and additional scalar fields S [47] . With 



the appropriate assignment of the horizontal quantum numbers to S, Higgs fields 4> Uj d and 
matter fields Q, d, u, one can construct non-renormalizable terms in the Lagrangian 



A 



1.1 



Q;O ll S'" -d i + -^g,o,,V" w/ ; + h.c. 



(3.1) 



which can give masses to the fermions if S has a vacuum expectation value (S). Here 
and riij are determined by the H charge assignment, so that the terms are invariant under H, 
M is a higher energy scale that communicates the horizontal symmetry breaking to the light 
states, and \fj U are some coefficients of order 1. The horizontal symmetry is spontaneously 
broken when S gets a vacuum expectation value, which introduces a small number e = (S) /M 
(this is the Froggatt-Nielson mechanism |55|]). Different powers of e in the Yukawa coupling 
account for mixing angles and the hierarchy of fermion masses. In supersymmetric theories, 
squarks have the same H charges as the quarks of the same multiplet and will obtain masses 
by the same mechanism. An astute choice of charges can allow for the alignment of squark 
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mass matrices with quark mass matrices, thereby suppressing flavor- changing effects. Notice 
that more than one U(l) symmetry is generally needed in order to get feasible models 
consistent with the experimental bounds. In p3| , the alignment model is associated with 
spontaneous CP-violation, which can predict small values of supersymmetric CP-violating 
phases so that the EDM bounds are satisfied. After diagonizing both the quark and squark 
mass matrices, we found that the (12) and (13) mixing angles in the gluino coupling vertices 
(although the (12) mixing angle is too small for this to be relevant) and the (13) component 
of the CKM matrix can have a CP- violating phase 0(1), while the (23) mixing is almost 
real, so there will be no sizable contribution to the dilepton charge asymmetry for B s . 

33, B7L B8L EH, 



The non-abelian models 135 



[1], f|2], [43] are motivated by the large 
top mass and the different behavior of the third family with respect to the first two families. 
The maximum flavor symmetry group is U(3) 5 in the absence of Yukawa couplings. It can 
be assumed that there is a non-abelian flavor symmetry Gf C U(3) 5 where the first two 
families and the third family transform differently; Gf can be continuous or discrete, gauged 
or global. There are a variety of models based on different Gf: Gf = U(2) |36], f40| , 
G f = SU(2) |3J, |TJ, Gf = (S 3 ) 3 @, 0, G f = U(l) x 0(2)/Z 2 @ and G f = A(75) @. 

and is broken by some small factor A. Because the 



This symmetry is only approximate, 
symmetry guarantees that the first two families are nearly degenerate, FCNC for the light 
generations are heavily suppressed. In some models where the degeneracy of the first two 
families does not fully resolve the FCNC constraints [5], [IT], [42], different scenarios are 



proposed to relax the constraints: the first two generations can be much heavier than the 
third generation m 12 ~ 10rh 3 (scenario (a) in [[5J), or the CP- violation phase is very small 
ip ~ 1(T 2 (scenario (b) in g(| and 

~T6f provide another possible solution to the FCNC and 



The heavy squark models 44, 45 



CP-violation problems by allowing the maximal masses consistent with naturalness bounds. 
All those models permit the first two generations of squarks to be heavier, which is crucial for 
solving FCNC problems. In effective supersymmetry [4T3|], a new gauge group G is introduced, 
which enlarges the accidental symmetry group and thus forbids the renormalizable B- and 
L-violating interactions. It also introduces a new mass scale M ~ 5-20 TeV, which sets the 
mass scale for the first two generations. In addition, the requirement of naturalness implies 
that some squarks and most gauginos must have a mass below ~ 1 TeV. In [j44| , 

the first two families are charged under a gauged U(l) symmetry, while the third family is 
neutral. Thus, there is an extra contribution to squark masses of the the first two families 
coming from the D-term, which generates the mass hierarchy. The constraints on the mixing 



angles come from the naturalness of the Higgs sector and the squark mass matrices [49 



We should emphasize that the heavy-squark models cannot fully satisfy the FCNC and CP 
constraints by themselves. They have to be combined with non-abelian symmetry (scenario 
(a) in [f4~0"| |) or have some alignment for the squark mass matrices. 

In the next section, we select specific models from the above-mentioned papers and 
present a detailed study of the dilepton asymmetry and how it can put constraints on the 
squark masses, mixing angles and CP-violation phases for these models. Of particular in- 
terest are those models that specified the squark and quark mass matrix textures (therefore 
we do not consider ||43||) and have large mixing with the third generation (so we do not con- 



10 



(23) mixing 



(13) mixing 



Model 



LL RR 



LL RR 



B 



|32 

A |3 

El 



]38 



B+C |0 



35 



A 2 A 4 LL > RR 
A 2 1 LL <C 
Small CP-violation angle 



m b 



"A 2 " 

A 2 
A 3 
A 2 



A 2 
A 5 
A 1 



LL < RR 
LL = RR 
LL > RR 
LL > RR 



A 3 A 3 LL = RR 

Too small mixing 
A 3 A 7 LL > 
"A^ 



A 3 
A 2 
A 3 



A 3 ^ 
A 3 
A 4 
A 4 



LL < L1R 
LL = L1R 
LL > RR 
LL ~ 



Mass 
insertion 



C p[ 



A 2 
"A 2 " 



A' 



LL = RR 



LL > 



A 3 
"A 7 " 



A a 



LL = RR 



LL > 



Vertex 
mixing 



Table 1: Selected models from the literature, which will be analysed in section |j. Here A 
alignment, B — non-abelian, C — heavy squarks. 



sider |39|]). We also require that the models we consider satisfy the FCNC and CP- violation 
bounds set by the experimental value of Ek, AMx, AMo, d e and d n . This excludes models 



in 



36|, fill ) an d bounds the squark masses in some other models. In 



m~ q > 200 GeV, 



so that the quark-squark alignment solution to Atuk will not run into problems with Amu. 
In the scenario (a) of p0| , ^q 3 ,d 3 should be heavier than 550 GeV if rriq ~ rrig, so that the 
super symmetric contribution to Ek is within the experimental bounds, but this constraint is 
relaxed if the CP- violation phase is small, so we do not impose this bound. In [5B|, there is 
1 TeV upper bound which comes from the naturalness of the Higgs sector. For the models 
where the squarks of the first two generations are heavy (scenario (a) in ^(| and j46]), m 12 
are required to be heavier than 10m 3 . 

The models should have an 0(1) CP-violation phase difference in the (13) and (23) 
mixing angles between the SM and supersymmetric models, so that their contributions to 
the dilepton asymmetry is not negligible. Thus, in [M\, we only consider (13) mixing since 
the (23) mixing CP-violation phase difference is small. 

In Table [l], we list different models with the predicted LL and RR mixing (up to 0(1) 
uncertainties). As will be shown in the next section, M 12 can get contributions from both 
LL and RR mixing. There are two different cases: 

1. Either LL or RR mixing dominates. 

2. LL mixing is comparable with RR mixing. For definiteness, we take LL = RR in our 

where the (13) LL and RR mixing are A 3 and 

RR in Table [l]. It is nonetheless important 



g2 



calculation, except for model in ref. 
A 4 respectively, which is denoted by LL 
to take the RR term into account because the contribution to M i2 from (#31) ll(^3i) rr 
is large. 

We also show in the table whether the mass insertion (with mixing parameter Sll^rr) or the 
vertex mixing (with mixing parameter K^ r) method is used. We only use the vertex mixing 
in scenario (a) of ^TJ and in , where the first two generations are much heavier than the 
third one, while the mass insertion should be a good approximation in the other cases, though 
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it should be noted that in a detailed analysis one would account for the nondegeneracy of 
the squarks. 



4 Numerical Results 



We will assume the measurement of Im(Ti2/Mi2) for either Bd or B s with a sensitivity 
2 x 1CT 3 . This measurement can be obtained from either the single or dilepton asymmetries. 
For defmiteness, since An ~ Im(r 12 /M 12 ), we refer to the dilepton asymmetry in this section. 

From Eq. (|2.26| ), we see that the dilepton asymmetry depends on (Ar/AM)sM, h and 
9. (Ar/AM) SM is given in Eq. fl2.27| ) which has been calculated in the SM; 9 is the phase 
difference between Mf 2 M and Mf 2 USY , which can be in the range from to 2n. Since \Au\ is 
symmetric with respect to n, we only consider 9 to be in the range of to tt. The quantity 
h is the ratio between the amplitudes of Mf 2 USY and Mf 2 M , which can be calculated through 
the AB = 2 box diagrams with q, g or q, W running in the loopQ; Mf 2 M has already been 
calculated, including QCD corrections 0, and gives (we only include the top-quark 

contribution since this is the largest effects): 



M 



SM 
12 



^ m 2 t B B JlM Bq (v tq v; b f A{zt] 



2f 



where 



1 9 
4 + IT 



z t 



2 (1 - z t ) 



3 z 2 \nz t 



VQCD, 



Z t 



m. 



m 



(4.1) 



(4.2) 



Bs q is the "bag" parameter describing the uncertainty in evaluation of the hadronic matrix 
element, Ms q and /#, are the B q meson mass and decay constant respectively. Although 
B Bq f B( MB q cancels in h, as will be shown below, we still need to know the value of BB d 
and fs d since we have to take into account the constraint from AM^. The recent values 
of B Bd and f Bd are B Bd = 1.29 ± 0.08 ± 0.06 [0, f Bd = 175 ± 25 MeV [[51]. t] QC b is the 



QCD correction factor, which is taken to be 0.55 ± 0.01 [2(| in our calculation. Mf 2 can 



be calculated either in the scenario of vertex mixing |3(| or using mass insertion p8 
first case, Mf 2 USY is given in terms of Kf R (we take (Kf R ) 33 ~ 1)0 



In the 



M 



VM 
12 



216m| 3 



l -B 



B q f< 



B a M Bq 



K 



d\2 
L)3i 



36-24 



M. 



B a 



2 - 2x + (1 + x)lnx 



+ (K d R )l i )(66U(x) + 2AxU(x)) + 



72 + 384 



M, 



B a 



(x-1) 



1 — x 2 + 2xlna; 



x 



(X 



1 I 1 



(4.3) 
(4.4) 



4 Here we do not consider the contribution from box diagrams with q and chargino (or neutralino) running 
in the loop since it is suppressed by a factor ((ft/ffs) 4 with respect to the q, g box diagram contribution. 



J There are mistakes in the formulas given by Ref. 
corrected formulas here. 



which is pointed out by Ref. |28|. We use the 
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f 6 (x) 



f 6 (1) = 1/20 




Figure 2: Plot of f±{x), /i(x), fe(x), fe,(x)with x in the range of 1/3-3. 



where i = 1, 2 for Bd and B s respectively. In the mass insertion notation, 
al 1 



mW = 

12 216m? 3 



-B Bq f B M Bq {{(ShYll + fe)i*)(66/ 6 (x) + 2Axh{x)) + 



RR 



36 - 24 



/ 6 (x)+ 72 + 384 



M Ba 



6(1 + 3a;)lnx + x 3 - 9x 2 - 9x + 17 
6(x- l) 5 ' 

6x(l + x)lnx — x 3 — 9x 2 + 9x + 1 
3(x - l) 5 ' 



The mass parameters we used in our calculation below is M Bd 
5369.3 ± 2.0 MeV, m w = 80.41 ± 0.10 GeV pi. From Eqs. 



(4.5) 
(4.6) 

(4.7) 

5279.2 ± 1.8 MeV, M Bs = 
ID, (O) and (Oh, we can 



see that /i depends only on rrig, mixing parameters, and x. The x dependence comes from 
the functions fi(x), f±{x) or fe(x), fe(x). We plot the x dependence of these functions 
in Fig. H] and also show their values when x = 1, which is the value of x we use in the 
calculation below. The supersymmetric contribution measured by h can dominate over the 
SM prediction for AM Bq because of large mixing angles or small m q . In Eq. (|2.26|) , we 
see that An is suppressed when h is either too large or too small; thus the limit on the 
value of An that will be experimentally accessible translates into a range of h that can be 
probed, given a specific 6. This in turn determines the range of rrig that can be probed in 
particular models, which specify (approximately) the mixing angles. Again we emphasize 
that in general for detectable values of 9 we will find both an upper and lower bound on h 
which will be experimentally accessible. The upper bound corresponds to too large mixing 
whereas the lower bound corresponds to too small a supersymmetric contribution (so too 
small CP violation). The upper limit on h translates into a lower bound on the squark mass 
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Vertex mixing Mass insertion 



6 LL (RR) dominates LL - RR LL (RR) dominates LL - RR 
tt/2 rrig -192-505 m 9 ~ -323-849 -45-119 m f -224-590 

3tt/4 m g - -192-485 m g - -323-814 m, -45-114 m g - -224-566 

Table 2: (GeV) range from dilepton asymmetry with An = 2 x 10~ 3 , x = m~/m~ = 1 
and the mixing angles are the same as the corresponding CKM entry. 

Vertex mixing Mass insertion 

6 LL (RR) dominates LL - RR LL (RR) dominates TT^RR 
tt/2 m 9 ~ -130-310 m~ q -214-509 m 9 ~ -31-73 m f -150-357 

3?r/4 mg -115-352 m g - -188-578 m g - -27-83 m q - -132-405 

Table 3: (GeV) range from 5 S dilepton asymmetry with An = 2 x 10~ 3 , a; = m~/m~ = 1 
and the mixing angles are the same as the corresponding CKM entry. 

(for definite mixing angles) whereas the lower bound translates into an upper bound on the 
mass. 



For our results we take into account the constraint AM^ = 0.470 ±0.019 ps [|16|]. The 
super symmetric contribution to AMs d cannot exceed this, which puts a severe bound on the 
(13) mixing and rrig. The B s mixing has no such constraint. 

As explained in Sec. ||, Bd and B s both contribute to the total dilepton rate. Assuming 
a measurement on the dilepton asymmetry with a sensitivity of 10~ 3 , we require at least 
2 x 10~ 3 dilepton asymmetry from either B^ or B s . This is also the current upper bound on 
the Bd standard model contribution, so it is clearly identifiable as new physics. Of course if 
better precision is possible a measurement of the B s system with greater accuracy would be 
interesting since its standard model rate is much lower. 

In Eq. ( p.28|) , we see that the phase difference between Mf^ SY and Mfg must be large 
enough so that the dilepton asymmetry can be measured. For the precision we choose, we 
find that 

f 30° -40° fori?,, 
~\ 52° - 106° for B s . K ' 

Here the large range in the lower limit of 9 comes from the large uncertainties in (Ar / AM)sm- 
Models with too small a CP-violation phase difference between the SM and supersymmetric 
models cannot be tested using the dilepton asymmetry. 



For the results presented below, we take the SM CKM matrix elements [16, 55 : 



\V td V* b \: 0.0066 - 0.0102, WMl ■ 0.026 - 0.060. (4.9) 

In Tables |2] and || we show the allowed ranges of rriq when the dilepton asymmetry of 
Bd or B s is larger than 2 x 10 -3 . Here we assume that the (13), (23) mixing angles are the 
same as the corresponding CKM entries. We notice that the vertex mixing results are larger 
than the mass insertion ones. This is because we only consider the contribution from the 
third generation in the vertex mixing case, while all the generations contribute in the mass 
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(23) mixing 



(13) mixing 



nig (GeV) 



mq (GeV) 



Model 




LL 


RR 


= tt/2 


a = 37T/4 


T T 

LL 


RR 


e = tt/2 


= 3tt/4 




|32 








\2 
A 


A 4 


20-112 


IS 1 07 


A 


A 3 


176-715 


176-686 


A 


33 


^ 


35 


a 


\2 
A 


1 


511-2807 


A Q1 S7 

4oU— olo ( 




Too small mixing 




34 










Small CP-violation 


angle 


A 3 


A 7 


36-144 


36-139 




m\ b 




A 2 


A 1/2 


229-1255 


201-1425 


A 3 


A 3 ^ 


397-1614 


397-1549 


B 


37 


■ 


40 1 b 


A 2 


A 2 


100-549 


88-624 


A 3 


A 3 


176-715 


175-686 




38 








A' 3 


A 5 


4-22 


4-25 


A 2 


A 1 


178-722 


178-693 




42 








A 2 


A 4 


20-112 


18-127 


A 3 


A 4 


73-298 


73-286 


B+C 


[40 


a 




A 2 


A 2 


143-784 


126-890 


A 3 


A 3 


253-1029 


253-987 


C 


46 








A 2 




87-476 


76-541 


A 3 




151-613 


151-588 



Table 4: nig ranges for different models with An = 2 x 10 3 , x = m~/ni~ = 1 and A = 0.2 
for B s ((23) mixing) and ((13) mixing) dilepton asymmetries. Notations are the same as 
in Table 0. 



insertion method. In the latter case, the leading-order contribution cancels because of the 
GIM mechanism, which gives smaller Mf 2 USY . Thus, the squark masses obtained using the 
mass insertion method have to be smaller in order to compensate for this weakening effect. 
However without a detailed knowledge of the mixing angles and masses either method must 
be viewed as an approximation. 

We need to take into account the current experimental lower bounds on the squark mass 
coming from the non-observation of any supersymmetry signals at either LEP |52] or the 
Tevatron [^, [54[]. For the vertex mixing case, when the first two generation squarks are much 
heavier than the third generation ones, the nig bounds are on the lightest sbottom mass. 
The sbottom masses are all larger than 110 GeV, which have not been excluded by the 
experimental lower limit (m^/2) p2]| . For the mass insertion case when the squark masses 
are almost degenerate, we impose the constraints from CDF and D0 searches. While LL (or 
RR) dominates, the squark masses are smaller than 120 GeV, which has already been ruled 
out by the CDF (nig > 230 GeV f53|) and the D0 (m 9 - > 260 GeV|4[) limit; this means that 



if the soft supersymmetry-breaking models have the squark mixing angles in the SM ranges 
and LL (or RR) is dominating, then the dilepton asymmetry is too small to be measured 
experimentally. If LL = RR, the squark masses can be in the range of 130-600 GeV; most 
of the mass ranges have not been ruled out yet. Thus, the corresponding supersymmetric 
models can still be tested through a dilepton charge asymmetry measurement. 

Also notice that the lower limits of nig in Table ^ (B^ case) are the same for different 9. 
This is because the upper bound on h (which corresponds to the lower bound on rriq) is set 
by the experimentally measured value of AM# d . 

We now consider the implications of our results. In Table |4], we show the ranges of nig 
for different soft supersymmetry-breaking models. The lower mass limits coming from the 
constraints of FCNC have been imposed. We see that the models divide into several groups: 
those that show an effect in both, neither, or one of the Bd and B s systems. It should of 
course be remembered that we used a definite value for the angle in these tables so there is an 
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order one fudge factor available in any model. Nonetheless, there are some clear tendencies 
indicated by these results. 

For the models in |37|, |4(J, both the (#13,23)1,1, and (#13, 23) rr are of the order of A 3 and 
A 2 , respectively, which are comparable to the SM CKM mixing angles. This is the optimal 
situation in that for reasonable supersymmetric masses, the standard and nonstandard box 
contributions are competitive. The squark masses are below 1 TeV, which can be explored 
experimentally. 



The models in |33], |35| predict quite different mixing angles from the SM. The RR mixing 
can be as large as 1 in (23) and 0(A 3 / 2 ) || in (13), which increases the squarks to heavy 



masses in order to keep the supersymmetric and SM box diagrams to be comparable to 
each other (h ~ 1). If the squark masses are indeed light, it is likely that there is a large 
supersymmetric contribution to flavor violation and that it will be completely undetected 
in the measurement of a lepton asymmetry! Since such large mixing angles will give rise to 
noticeable effects in other measurements, the lepton asymmetry provides a nice complement 
to such measurements and would provide a clear determination of a model with large mixing 
angles. However, it is possible that the squark masses are at the lighter end (which would 
be measured) and that there is visible evidence of a lepton asymmetry. 

There are several models which should give a reasonable asymmetry for Bd but not for 



B s . These include the models in |32]. |38], 4"2"| , although the latter only has a small range 



of squark mass which has not already been excluded by the CDF and D0 limits. Actually 
for pij , the squark mass range which is accessible is already at the limit of what would be 
permitted once the 200 GeV lower bound arising from the Amp constraint is accounted for. 
We can see that the model of Ref. P3 which explicitly addresses the CP violation as well 



as the flavor structure gives no measurable effects in the lepton asymmetries. 

We also find that the heavy squark models are very likely testable. However, we note that 
we assumed the mixing angles agreed with their CKM counterparts. Should they be bigger, 
the asymmetry could be too small to measure. For the effective supersymmetry model in 
H| , the upper limit increases with the mixing angles; it is 600 GeV for the mixing to be in 



the SM ranges. If the mixing is too large, this model predicts too small dilepton asymmetry 
because of the 1 TeV upper limit in squark masses coming from the naturalness constraints 
and thus cannot be tested. 

In the analysis above, we assumed that the dilepton asymmetry of 2 x 10 -3 is measured for 
either Bd or B s . If we assume 4 x 10~ 3 sensitivity, then the lower limit on 9 is 56° — 72° for Bd 
and 89° — 139° for B s . For any given experimental sensitivity D, we define a ^-dependent 
scale factor Sd, s (@, 2x10^ )' ^ ne sc L ua rk mass under this new sensitivity is then (rrig) D = 
Sd,s{&, 2xio~ s )( m g)2xio~ 3 - If we take $ = 37r/4 and D = 4 x 10 -3 as an example, in Bd, 
the lower limit remains the same because it comes from the AM# d constraint, while for the 
upper limit 5*^(3^/4, 2) upper = 0.765. In B s , the reduction in squark mass range follows from 
5 s (3?r/4, 2) U p per = 0.78, S s (3n/4:, 2) down = 1.28. In general, the factor D can be obtained 
from Eq. ( gg| ). 

We conclude that models are certainly distinguishable from the lepton asymmetry mea- 
surements alone. With other complementary measures there is some hope to resolve the 
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flavor problem of supersymmetry. However, we have assumed a reasonably good sensitivity, 
which is essential for the measurement to be useful. 



5 Conclusions 



It is possible that the dilepton asymmetry could be one of the first indications of physics 
beyond the standard model. Once the source of new physics is ascertained through direct 
measurements, it can be used to impose further constraints. In particular, if the new physics 
is indeed low-energy supersymmetry, the dilepton charge asymmetry can be used to distin- 
guish various soft supersymmetry-breaking models. The range of parameters which can be 
tested is such that there is often a good overlap with interesting flavor models. Unfortu- 
nately, there is not an unambiguous identification of the size of the signal with the category 
of model; nonetheless particular models with definite patterns for masses and mixings can be 
tested. We emphasize the importance of taking a model-dependent approach; although large 
mixing angles can give big effects in searches for physics beyond the standard model, existing 
constraints and the attempt to motivate flavor physics parameters through an underlying 
symmetry structure often disfavor this assumption. 

We emphasize that there are uncertainties in the precise angles which can change the exact 
range of squark mass which is covered. In particular, all models we presented have order unity 
uncertainties for the angles which can change the range of mass which is probed. Furthermore 
nondegeneracy of the squark mass and gluino mass introduces another parameter which can 
affect the precise range of parameters which is covered. Finally, the mass insertion method 
is an approximation; it is generally only a good one if the mass parameter is appropriately 
interpreted. Nonetheless, the overall message is clear; it would be very interesting to do an 
accurate measurement of the single lepton and/or dilepton asymmetries. 

In this paper, we have focussed on the lepton charge asymmetry as a means of searching 
for new physics, in particular for testing new models of flavor. These lepton asymmetries 
are sensitive to new sources of mixing in the B system. These new sources of mixing can 
be independently tested in other measurements |], [TJ|, [1£], [56], pTFfl . Particularly interesting 
alternative tests for CP- violation in B s are the study of B s — > DfD~, B s — > tp<f), and 
related modes JT7], E8|, |51| ; B d , in particular modes such as B d — > (pK s and B d — > K s tt° 
|57| , [59], where new physics effects might be large, and measurements of deviations in 
precisely predicted rates such as B — > J/ipK s . In particular, models with large mixing in 
the B s sector (of order unity) should give large effects for these measurments. Even with 
these studies, it will be useful to have an alternative means of searching for new physics. 
Precisely because the rate is negligible in the standard model, the lepton asymmetry would 
be an extremely important measurement. 

Acknowledgements: We are grateful to U. Aglietti, M. Lusignoli, G. Martinelli, A. Nelson, 
and especially Y. Nir, R. Rattazzi, D. Wood, and H. Yamamoto for useful conversations. 
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